Thermodynamic properties of the coupled dimer system Cu 2 (C5Hi 2 N 2 )2Cl4 



o 
o 



S. Capponi and D. Poilblanc 

Laboratoire de Physique Theorique, Universite Paul Sabatier, CNRS, 31062 Toulouse, France. 

(Dated: February 6, 2008) 

We re-examine the thermodynamic properties of the coupled dimer system Cu2(C5Hi2N2)2Cl4 
under magnetic field in the light of recent NMR experiments [Clemancey et al., Phys. Rev. Lett. 
97, f 67204 (2006)] suggesting the existence of a finite Dzyaloshinskii-Moriya interaction. We show 
that including such a spin anisotropy greatly improves the fit of the magnetization curve and gives 
the correct trend of the insofar unexplained anomalous behavior of the specific heat in magnetic 
field at low temperature. 

PACS numbers: 75.40.Cx,75.10. Jm 
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The molecular solid Cu2(C5Hi2N2)2Ci4 is of particular 
interest due to the formation of dimers between the two 
copper (S=l/2) spins of each molecule. It was long be- 
lieved that the coupling between those dimers were quasi- 
one dimensional hence realizing the physics of a two-leg 
spin ladder jlj. Later, it was emphasized that the mag- 
netic coupling between molecules might be more three- 
dimensional (3D) 0]. In a recent Letter, Clemancey 
et al. 3] proposed that including Dzyaloshinskii-Moriya 
(DM) interactions is essential to explain its experimental 
behavior under magnetic field. From simple symmetry 
considerations, the lack of inversion centers at the mid- 
dle of the dimers enables a priori the existence of such 
a DM term which can lead to experimentally observable 
effects such as the appearance of a uniform transverse 
magnetization [4[ as shown by T=0 numerical calcula- 
tions. However, a theoretical investigation of physical 
properties at finite T in the presence of such a DM inter- 
action has not been done sofar. 

Following 0] we consider an anisotropic spin ladder 
under a magnetic field h along z with a staggered DM 
term along y, 

T~i = J± J]] • S^2 + J\\ / ^(Si,i • S i+11 + Si 2 • S i+12 ) 

i i 

+ ^(-iy y .(s. M x s ii2 ) + s|lB ^(^ + sy, 

i i 

where the last term corresponds to the field Zeeman en- 
ergy. Our aim here is to show that the thermodynamic 
properties of the above-mentioned molecular system can 
be partly explained within this model without invoking 
more complicated three-dimensional couplings. 

Specific heat measurements done on the same com- 
pound could not be satisfactorily explained with a SU(2)- 
invariant Heisenberg interaction i.e. D = d, H, 0]- 
Indeed, despite a fairly good agreement at low fields, a 
large discrepancy occurs for fields around 10 T for which 
a strong anomalous low-temperature peak around 1 K 
is seen experimentally and not reproduced theoretically. 
In contrast, we show here that the inclusion of the fi- 
nite DM interaction of Q provides a more quantitative 
agreement with the experimental data. We take param- 



eters of previous literature, J± = 13.1 K, Jm = Jj_/5 and 
g = 2.1, and use Exact Diagonalisations (ED) to compute 
the magnetization curve and the specific heat C mag (T, h) . 
Interestingly, a finite D strongly reduces finite-size effects 
so that accurate estimates can be obtained on 16-site sys- 
tems (see later). 
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FIG. 1: (Color online) Uniform magnetization per site m z 
vs magnetic field h (in Tesla) for various cluster sizes, (a) 
D = : m 2 exhibits finite-size plateaus at rational values 
multiple of 1/N that disappear in the thermodynamic limit, 
(b) Finite D : we obtain continuous curves (since S z is not 
conserved any more) which show very small finite-size effects. 
The square root singularity disappears and the magnetization 
curves become smooth. 

We start by computing the zero-temperature uniform 
magnetization along the field m z for different sizes. As 
is shown on Fig. [TJa), when D = 0, we recover the typ- 
ical behavior expected for an isotropic interaction with 
square root singularities |8| at the two critical fields hci 
and hc2- I n that case, since S z is conserved, it is suffi- 
cient to compute the ground-state energy in all S z sectors 
(in the absence of any magnetic field) and then perform 
a Legendre transform to get m z . As a side remark, we 
can mimic the thermodynamic limit by drawing a line 
connecting the middles of the finite-size plateaus of the 
largest available system (N — 32 sites). In contrast, for 
any finite D (Fig. O^b)), since S z is no longer a good 
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FIG. 2: (Color online) Normalized magnetization along mag- 
netic field m z : experimental data are taken from Ref . |5J nu- 
merical results are presented for various D (in units of J±) 
and N = 20 spins (lines are guides to the eye). In the specific 
case D = 0, we plot an extrapolation (see text). The best 
agreement with experiment is obtained for D ~ 0.05Jj_. 



quantum number, the magnetization has to be computed 
for each h value, which makes the problem harder. For- 
tunately, as shown on Fig. [ljb), the absence of SU(2) 
symmetry makes the finite-size effects almost negligible 
for any finite D. Moreover, we observe that the absence 
of singularities at hex and hc2 is well reproduced by any 
finite D as noticed in Ref. 0. 

Prior to the calculation of C(T, h), we estimate D from 
a fit of the experimental magnetization curve [9( at low 
temperature. On Fig. [2J we plot the numerical uniform 
magnetization for various D. By comparing to experi- 
mental data, we find that D ~ 0.05 J± gives an excellent 
fit with no other adjustable parameter. Remarkably, our 
estimate of D is exactly the numerical value obtained in 
Ref. S 

We show on Fig.[3]the main results of our specific heat 
calculations and compare it to the experimental data of 
Ref. @. First, at low magnetic field up to 6 T (h < h c i), 
D has almost no effect due to the very large spin gap 
(the D = and D = 0.05 J± curves are indistinguishable 
in Fig(3]Ja)) and the agreement is very good. However, 
above this magnetic field, we observe, in addition to the 
large broad maximum around 6 K, another peak around 1 
or 2 K. In our results, this low-temperature peak has the 
largest magnitude around h ~ 10-11 T, as seen in exper- 
iment. By taking our previous value, D/J± = 0.05, we 
clearly observe a significant increase of its intensity that 
translates into a better comparison with experiments. 

It is of interest to see whether a larger D value could 
give rise to a better agreement for the heat capacity (at 
the price of worsening the magnetization fit). However, 
as indicated on Fig. [H we observe that increasing the 
DM interaction has two important effects : first, it leads 
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FIG. 3: (Color online) Magnetic specific heat vs temperature 
for various magnetic fields h (in Tesla) : experimental data 
(symbols) are taken from Ref. @; exact numerical results are 
presented for D = (dashed lines, obtained with 20 spins) 
and D — 0.05Jx (full lines, obtained with 16 spins) showing 
a better agreement. 



to an increase of the height of the low-temperature peak 
but, secondly, to a shift of its position in temperature 
with little influence on the second (broad) maximum. As 
a consequence, the minimum between the two peaks in- 
creases too with D. 

Therefore, a small DM term can partly account for the 
anomalous experimental height of the low-temperature 
peak. Note that small sources of discrepancy between 
theory and experiment should still remain like the exper- 
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imental uncertainty (due to the subtraction of the phonon 
contribution) , the deviation of the compound from a real 
ladder system, the small finite size effects in our calcula- 
tion, etc... 
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FIG. 4: (Color online) Magnetic specific heat for h — 10 T 
obtained numerically for various D/J± = 0, 0.05, 0.1 and 
0.15 with 16 spins. As D increases, we observe that both the 
low-temperature peak intensity and position increase. 



Lastly, we would like to discuss briefly the finite size 
effects at finite T on e.g. C mag (T) by comparing data 
obtained on clusters with 12 and 16 spins. As shown in 
Fig. 03a), in the absence of a magnetic field, finite size 
effects are invisible on the scale of the plot. Moreover, the 
specific heat is almost unsensitive to a small finite DM 
term. This contrasts to the case of intermediate h where, 
as shown above on Fig. 2J the data are greatly sensitive 
to the value of D. In addition, Fig. E^b) clearly shows 
that, fortunately, finite size effects are reduced when D = 
0.05Ji (compared to D = 0), the two curves for N = 12 
and N = 16 being almost superposed. This provides a 
great deal of confidence on our previous analysis. 



To conclude, we have shown that the proposal of Ref.0 
of a finite (small) DM interaction gives the correct trend 
of the insofar unexplained anomalous behavior of the spe- 
cific heat in magnetic field of Cu2(C5Hi2N2)2Cl4. From 
magnetization measurements, our considerations provide 
an estimate of D around 0.05Jl in agreement with the 
value extracted from other recent measurements 
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FIG. 5: (Color online) Magnetic specific heat for h = (a) 
and h = 9 T (b) obtained numerically for D/J±_ = (full and 
dashed lines) and 0.05 (dot-dashed lines) with N = 12 and 
N = 16 spins. Note that in (a) all four curves are indistin- 
guishable. 
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